THE K-MOMENT PROBLEM FOR CONTINUOUS LINEAR 

FUNCTION ALS 



JEAN B. LASSERRE 

Abstract. Given a closed (and non necessarily compact) basic semi- 
algebraic set K C R n , we solve the K-moment problem for continuous 
linear functionals. Namely, we introduce a weighted £i-norm £ w on 
R[x], and show that the ^ w -closures of the preordering P and quadratic 
module Q (associated with the generators of K) is the cone Psd(K) 
of polynomials nonnegative on K. We also prove that P and Q solve 
the K-moment problem for £ w -continuous linear functionals and com- 
pletely characterize those f w -continuous linear functionals nonnegative 
on P and Q (hence on Psd(K)). When K has a nonempty interior we 
also provide in explicit form a canonical ^-projection gj for any poly- 
nomial /, on the (degree-truncated) preordering or quadratic module. 
Remarkably, the support of gj is very sparse and does not depend on 
K! This enables us to provide an explicit Positivstellensatz on K. At 
last but not least, we provide a simple characterization of polynomials 
nonnegative on K, which is crucial in proving the above results. 



1. Introduction 

This paper is concerned with basic closed semi-algebraic sets K C M. n 
and the preordering P(g) and quadratic module Q(g) associated with the 
finite family of polynomials (gj), j £ J, that generate K. In particular, 
when K = W 1 then the latter two convex cones coincide with the cone 
S[x] of sums of squares (s.o.s.) of polynomials. The convex cones P(g) 
and Q(g) (which are subcones of the convex cone Psd(K) of polynomials 
nonnegative on K) are of practical importance because on the one hand 
nonnegative polynomials are ubiquitous but on the other hand, polynomials 
in P(g) or Q(g) are much easier to handle. For instance, and in contrast 
with nonnegative polynomials, checking whether a given polynomial is in 
Pd{g) (C P(g)) or Qd{g) (C Q{g)) (i.e., with an a priori degree bound d on 
its representation) can be done efficiently by solving a semidefinite program, 
a powerful technique of convex optimization. 

The celebrated K-moment problem is concerned with characterizing all 
real sequences y = (y a ), a £ N ra , which can be realized as the moment 
sequence of some finite Borel measure on K. For such sequences y, the 
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linear form L y associated with y is called a ~K-moment functional and Hav- 
iland's theorem states that L y is a K-moment functional if and only if 
L y (f) > for all / G Psd(K). And so one says that P(g) (resp. Q(g)) 
solves the K-moment problem if the K-moment functionals are those which 
satisfy L y (f) > on P(g) (resp. Q(g))- However, this is true if and 
only if P(g) (resp. Q(g)) is dense in Psd(K) for the finest locally con- 
vex topology on M[x]. When K is compact, the K-moment problem was 
completely solved (with P(g) = Psd(K)) in Schmiidgen [15j and soon re- 
fined to Q(g) = Psd(K) in Putinar |13j for Archimedean quadratic modules 

Q(g). 

Since then, recent contributions have tried to better understand (in even 
a more general framework) the links between Psd(K) and closures (and 
sequential closures) of preorderings and quadratic modules, one important 
goal being to extend (or provide analogues of) Schmiidgen and Putinar's 
Positivstellensatze [15 1. \13\ to cases where K is not compact. In particular, 
Scheiderer |14j has shown rather negative results in this direction. For more 
details on those recent results, the interested is referred to e.g. Powers 
and Scheiderer |12| . Kuhlmann and Marshall [8], Kuhlmann et al. [9], and 
Cimpric et al. [5]. 

On the one hand, all linear functionals are continuous in the finest locally 
convex topology, on the other hand, the negative results of Scheiderer [14] 
suggest that solving the K-moment problem via preorderings or quadratic 
modules is possible only in specific cases, and so this topology is not the most 
appropriate in general. So why not rather consider other topologies on R[x] 
and the K-moment problem for linear functionals L y that are continuous 
for these topologies? This is the point of view taken in Ghasemi et al. [6] 
where the authors consider certain (weighted) norm-topologies and show 
that the closure of the cone of sums of squares is Psd([— 1, 1]™) as did Berg 
[3] for the £i-norm. Notice that this was also the point of view taken in 
Schmiidgen |16| for a class of non commutative (enveloping) *-algebra; in 
the latter context [16], the author proves that the cone of strongly positive 
elements is the closure of the (smallest) cone of sums of squares in certain 
topologies. 

Contribution. In view of the negative results in Scheiderer [14], we 
also consider the above mentioned viewpoint and look at the K-moment 
problem by using a particular weighted ^i-norm on R[x] (denoted £ w for a 
certain sequence w : N n —> M>o) rather than the usual finest locally convex 
topology. In this framework we solve the K-moment problem for basic closed 
semi-algebraic sets K C W 1 in the following sense. We prove that (a) the 
^■-closure of P(g) and Q{g) is exactly Psd(K), and (b) P(g) and Q(g) solve 
the moment problem, i.e., the K-moment (£ w -continuous) functionals are 
those L y that are ^-continuous and nonnegative on Q(g) (or P(g)). In 
fact, such linear functionals L y are characterized by: 

- Ly(h 2 gj) > for all h G M[x], and every generator gj of K. 
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- 3M > such that \y a \ < Mw a for all a G N n . 

• Next, when K has a nonempty interior, there exist £ w -projections of a 
polynomial onto Pd(g) and Qd(g), where Pd(g) (resp. Qd(g)) denotes the 
subcone of elements of P(g) (resp. of Q(g)) which have a degree bound d 
in their representation. In general these projections are not unique but we 
provide a canonical £ w -projection gj for any polynomial /, which takes a 
remarkably simple and "sparse" form, and particularly when £ w is the usual 
^i-norm, in which case 

n 

(1-1) 9f = f + \ + J2 X i x i d > 

i=l 

for some nonnegative vector A G R n+1 . In other words, the support ||<7/||o of 
gj does not depend on K and does not depend on d either! The dependence 
of gj on the </j's that define K is only through the coefficients (A*|). In 
addition, the support is very sparse since ||g/||o < ||/||o+n+l- This confirms 
the property of the ^i-norm with respect to sparsity, already observed in 
other contexts. Minimizing the ^i-norm aims at finding a solution with 
small support (where ||x||o = #{i : Xi ^ 0}). Finally, the vector A in (jl.ip 
is an optimal solution of an explicit semidefinite program, and so can be 
computed efficiently. 

• We also provide a canonical ^-projection of / onto P{g) fl M[x]2d which 
is again of the form , and use this result to characterize the sequential 
closure P(g)^ of P(g) for the finest locally convex topology. Namely, 

P{gf = |/GR[x] : 3d s.t. Ve > 0, / + e^l + J>fj G P{g) J , 

and the same statement is true for the quadratic module Q(g). This latter 
result exhibits the particularly simple form q := (1 + Y^l=i x i) f° r ^ ne 
possible polynomials q in the characterization of P(g)^ (and Q(g)^) provided 
in e.g. [5j [9]; e.g., in [9] it is stated that one may take the polynomial 
(1 + |[x|| 2 ) s for some s. 

• Thanks to the characterization of canonical £ w -projections, we finally 
obtain a Positivstellensatz on K of the following form: / G Psd(K) if and 
only if for every e > 0, there is some d G N such that the polynomial 

/ + €(1 + E?=i Etl *f /(2*00 ^ in P(g) (or Q(g)). 

• At last but not least, and crucial for the above results, we prove a 
result of independent interest, concerned with sequences y = (y a ), a G N n , 
that have a finite representing Borel measure fi. Namely, we prove that a 
polynomial / is nonnegative on the support of /j, if and only if J h 2 fdfi > 
for all h GK[x]. 

The paper is organized as follows. In Section 2, and after introducing 
the notation and definitions, we present the intermediate result mentioned 
above. In Section 3 we show that P{g) and Q{g) solve the K-moment prob- 
lem for £ w -continuous linear functionals. In section 4, we provide explicit 
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expressions for the canonical £ w -projections onto P(g), Q(g) and their trun- 
cated versions. Moreover, we characterize the sequential closures Q(g)$ and 
P(g)* and we end up with a Positivstellensatz for K. 

2. Notation, definitions and preliminaries 

2.1. Notation and definitions. The notation B stands for the Borel a- 
field of W l . Let M[x] (resp. M[x]^) denote the ring of real polynomials in 
the variables x = (x\, . . . ,x n ) (resp. polynomials of degree at most d), 
whereas S[x] (resp. S[x]^) denotes its subset of sums of squares (s.o.s.) 
polynomials (resp. of s.o.s. of degree at most 2d). For every a G N n the 
notation x a stands for the monomial x" 1 ■ ■ ■ x" n , \a\ stands for the integer 

a := cci H \- a n , and a\ stands for the integer a! For an arbitrary set S C 

M. n , let Psd(S) denote the convex cone of polynomials that are nonnegative 
on S. 

For every i G N, let := {/3 G N n : Y^jPj ^ d } whose cardinal is 
s(d) = ( n ^)- A polynomial / G R[x] is written 

x^/(x) = /« x °' 

and / can be identified with its vector of coefficients f = (f a ) in the canonical 
basis (x Q ), a G N n . The support of / G R[x] is the set {a G N n : / Q + 0} 
and let ||/|| := card {a : f a ^ 0}. Denote by ||/||i the £i-norm J2 a \ f a \ of 
the coefficient vector f . 

Crucial in the sequel is the use of the following ^ w -norm which is a 
weighted £i-norm defined from the sequence w = (w a ), a G N n , where 
w a := (2|"|a|/2~|)!. Namely, denote by ||/|| w the £ w -norm ^2 a w a \f a \ of the 
coefficient vector f; hence the £i-norm corresponds to the case where w a = 1 
for all a G N n . Both t\ and £ w also define a norm on R[x]d. 

Let S p C W xp denote the space of real p x p symmetric matrices. For 
any two matrices A,B G S p , the notation A y (resp. >- 0) stands for 
A is positive semidefinite (resp. positive definite), and the notation (A,B) 
stands for trace AB. 

Let v d (x) = (x a ), a G NJj, and let B° G R s ( d ) xs ( d ) be real symmetric 
matrices such that 

(2.1) v d (x)v d (x) T = £ X ° B - 

Recall that a polynomial g G Mfx^d is a s.o.s. if and only if there exists a 
real positive semidefinite matrix X G R s ( d ) xs ( d ) such that 

g a = (X,B° a ), VaeN^. 

Let gj G M[x], j = 0,1,..., m, with go being the constant polynomial 
<7o( x ) = 1 f° r an x ) an d let K C 1" be the basic closed semi algebraic set: 

(2.2) K := {x G R n : #(x) > 0, j = 1, . . . , m}, 



THE K-MOMENT PROBLEM FOR CONTINUOUS LINEAR FUNCTION ALS 



For every J C {1, . . . , m} let gj := Ylkejdki with the convention g$ : = 1, 
and let vj := [(deggj)/2) (with vj := v^y). 

Definition 2.1. With K as in [EE). letP(g),Q(g) C M[x] and P k (g),Q k (g) C 
M[x]2fc be the convex cones: 



P(g) 
P k {g) 

Q{g) 
Qk{g) 




o-jgj : o-j G £[x], J C {1, . . . ,m} 



^ vjgj : o-j £ T,[x] k _ vj , 

k JC{1 m} 



JC{l,...,m} 



5^<T j5j - : crj G £[x], j = 1, 



i=o 



ffj ffj : G S[x]fc_ V;j , j = 1, . . . , 

3=0 



m 



The set P(g) (resp. Q(g)) is a convex cone called the preordering (resp. 
the quadratic module) associated with the gj's. Obviously, if h G P(g) 
(resp. h G Q{g)), the associated s.o.s. weights afs (resp. c/s) of its 
representation provide a certificate of nonnegativity of h on K. The convex 
cone P k (g) (resp. Q k {g)) is the subset of elements h G (resp. /i G Q{g)) 
with a degree bound 2k certificate. Observe that P k (g) C -P(<?) nRjxjjj. and 

Q fc ( ff )cQ( 9 )ni[x] 2fc . 

Moment matrix. With a sequence y = (y a ), ct G N ra , let L y : R[x] — > M 
be the linear functional 

/ (=^/ a x Q ) m- L y (/) = ^/ a y Q , /£%]. 

With d G N, the <i-moment matrix associated with y is the real symmetric 
matrix M^(y) with rows and columns indexed in NR, and defined by: 

(2.3) M d (y)(a,(3) := L y (x Q +^) = y a+[S , Va,/3GN2. 

Alternatively, M^(y) = X] a eN™ d 2/oBq. It is straightforward to check that 

{L y (g 2 )>0 V S GM[x] d } o M d (y) >: 0, d = 0, 1, . . . . 

A sequence y = (y a ) has a representing measure if there exists a finite Borel 
measure fi on M n , such that y a = J x. a dfi for every a G N n . Moreover, 
the measure \jl is said to be determinate if it is the unique such measure. 
Notice that with the £ w -norm on M[x] is associated a dual norm || • ||^ 
on the dual space M[x]* of ^-continuous linear functionals on R[x], by 
H-Lyllw = sup{|y a |/w a : a G N n }. 
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Localizing matrix. With y as above, J C {1, . . . , m}, and gj G R[x] (with 
gj(x) = ^ 7 gj-yx" 1 ), the localizing matrix of order d associated with y and 
gj is the real symmetric matrix Md(gjy) with rows and columns indexed 
by N^, and whose entry (a, f3) is just 
(2.4) 

M d (y)(gjy)(a,f3) := L y ( 5j (x)x a+ ^) = ^ 5 j 7 y a+/3+7 , Va,/3 G NJ. 

7 

If e S S W is defined by: 

(2.5) 9j (x)v d (x)v d (x) T = ^ B^x Q , Vxel", 

then M d (pjy) = ^ aG N™ Alternatively M d (gjy) = M d (z) 

where z = (z Q ), a € N ra , with z Q = L y (gjx a ). 

Multivariate Carleman's condition. Let y = (y a ), a G N n , be such 
that M d (y) >z for all del If for every i = 1, . . . , n, 

oo 

(2.6) £ L y(*fr 1/2 * = oo, 

fc=i 

then y has a finite representing Borel measure /U on M n , which in addition, 
is determinate; see e.g. Berg [3]. 

Closures. For a set A C R[x] we denote by A the closure of A for the finest 
locally convex topology on R[x] (treated as a real vector space). With this 
topology, every finite-dimensional subspace of R[x] inherits the euclidean 
topology, so that A also denotes the usual euclidean closure of a subset 
A C R[x](j. Following Cimpric et al. [5 J and Kuhlmann et al. [9], we also 
denote by A^ the set of all elements of R[x] which are expressible as the 
limit of some sequence of elements of A, and so A$ is called the sequential 
closure of A, and clearly A C A^ C A. If A C R[x] is a convex cone 

A X = {/ G R[x] : 3? G R[x] s.t. / + e q G A, Ve > }, 

and in fact, q can be chosen to be in A. Moreover, if A has nonempty interior 
(equivalently, has an algebraic interior) then A$ = A. 

Semidefinite programming. A semidefinite program is a convex (more 
precisely convex conic) optimization problem of the form 

inf { (C,X) : AX = b; X h }, 

for some real symmetric matrices C, X G S p , vector b G R m , and some linear 
mapping A : S p — > R m . Semidefinite programming is a powerful technique 
of convex optimization, ubiquitous in many areas. A semidefinite program 
can be solved efficiently and even in time polynomial in the input size of 
the problem, for fixed arbitrary precision. For more details the interested 
reader is referred to e.g. [17] . 
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2.2. A preliminary result of independent interest. Recall that in a 
complete separable metric space X, the support of a finite Borel measure 
\i (denoted supp \i) is the unique smallest closed set K C X such that 
MX\K) = 0. 

Theorem 2.2. Let f £ M[x] and fi be a finite Borel measure with all mo- 
ments y = (y a ), ol £ N n , finite and such that for some M > and all k € N 
and alli = l,...,n, L y (x 2k ) < (2k)\ M. Then: 

(2.7) />0 on supp fi ^> Jh 2 fdfi>0 V7t € R[x] 

M d (fy) bO, Vd = 0,l,... 

Proof. The implication => is clear. For the reverse implication, consider 
the signed Borel measure v{B) := j B fdfj,, for all Borel sets B € B, and let 
z = (z a ), a € N n , be its sequence of moments. By Lemma [5.1l the sequence z 
satisfies Carleman's condition (12. 6h . Next, recalling that Mfc(/y) = M^(z) 
for every k E N, 

Jh 2 fdfi>0 V/iGR[x]J ^ M fc (z) b 0,VfceN. 

This combined with the fact that z satisfies Carleman's condition yields 
that z is the moment sequence of a finite Borel measure tp on K n , which in 
addition, is determinate. Therefore, 

(2.8) z a = J x a /(x)d/i(x) = y x a dV(x), VaeN". 

cij/(x) 

Let T+ := {x : /(x) > 0}, r - := {x : /(x) < 0} and let = fi + + fi~ 
with /j, + (B) = (i(B n r+), H~{B) = n(B flT), for all B £ B. Similarly, 
let = i/^" — v~ with and ^~ being the positive measures defined by 
v + {B) = J B fdfi + and u~(B) = — J B fd/i~, for all B £ B. Since < 
fj,, one has 

J x 2k d[i + (x) < J x 2k d^(x) and / x} k dpC(x) < J x 2k dn(x), 

for all z = 1, . . . , n and all fcgN. Therefore, again by Lemma 15.11 both u + 
and v~ satisfy Carleman's condition (|2.6p so that both are determinate. On 
the other hand, (|2.8p can be rewritten, 

/x«^,x)=/x« rf ,-,x) + /x«*,x). V a£ N», 

and so = v + ip because v + and v + ip are determinate. But then 
= z^ + (r _ ) > z^~(r~) implies that v~ = 0, i.e., v = v + = t/j, and so the 
signed Borel measure v is in fact a positive measure. This in turn implies 
that / > for all x G supp// \ G, where G C supp/i and /z(G) = 0. 
Notice that by minimality of the support , supp fi\G = supp \i. Hence let 
x 6 supp// be fixed, arbitrary. As supp/i\ G = supp//, there is sequence 
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(xf) C supp fi\G such that — > x as I — > oo, and /(x^) > for all £. But 
then continuity of / yields that /(x) > 0. □ 

Interestingly, as we next see, Theorem 12.21 yields alternative characteriza- 
tions of the cone Psd(S) for an arbitrary closed set S C 1™. 

For a finite Borel measure \x (with all moments finite) and a polynomial 
/ G R[x], let [if be the finite signed Borel measure defined by ^f(B) := 
j B fd/jL for all B G B. Let 8^ := {// CT : a G E[x]}, i.e., 0^ is the set of finite 
Borel measures absolutely continuous with respect to fj,, and whose density 
(or Radon Nikodym derivative) is a sum of squares of polynomials. 

Let S[x]* C R[x]* be the dual cone of the cone of S[x], i.e., the set of 
linear forms L y on R[x] such that L y (h 2 ) > for all h G R[x], and similarly, 
let 9* C R[x] be the dual cone of &„, i.e., 6* := {h G R[x] : / ' hdv > 
0, Vz/ G (-)„}, 

Corollary 2.3. Lei S C R n 6e an arbitrary closed set and let [i be any finite 
Borel measure such that supp/i = S and x (-)• exp(|xj|) is [i-integrable for 
all i = 1, . . . , n. XTien wii/i / G R[x], 

(2.9) / G Psd(S) fif G S[x]* 

(2.10) Psd(S) = 9*. 

Proof. Let y = (y a ), a G N n , be the moment sequence of fj,. Observe that 
for every i = 1, . . . , n, and all k G N, 

L y((fb) ^ ^exp(|x<|)d/i(x) < M, 

for some M > 0. Moreover, supp/U = S and so by Theorem 12. 2| / > 
on S if and only if J s h 2 fd/j, > for all h G R[x]. Equivalently, if and 
only if f s o~d[if > for all a G S[x] (which yields (|2.9|) ). or if and only if 
Js /d/V > for all a G £[x], which yields (l2lH . □ 

3. The K-moment problem for ^-continuous linear functionals 

We first show that Q(g) (and P(g)) solve the K-moment problem for 
£ w -continuous linear functionals. 

When equipped with the £ w -norm, we may and will identify R[x] as the 
subspace of sequences with finite support, in the Banach space of real in- 
finite sequences f = (f a ), a G N ra , that are w-summable, i.e., such that 

T,a W a\fa\ < + 0O. 

Proposition 3.1. The dual of (R[x], || • || w ) is the space (R[x]*, || • ||^) of 
linear functionals L y associated with the sequence y = (y a ), a G N ra , which 
satisfy ||L y ||^ < oo, where ||L y ||w := sup{\y Q \/w a : a G N n }. 

Proof. If L y G R[x]* satisfies ||L y ||^ < oo, then 

\L y (f)\ < ^IfalWa^- < ||/||w||L y ||;, 
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and so L y is bounded, hence £ w -continuous. Conversely, if L y is ^-continuous, 
then consider the sequence of polynomials (f a ) C K[x] with x i— > / Q (x) : = 
x a /w a for every a G N n . Then ||/ a || w = 1 for every a G N n , and L y (f a ) = 
Ua/wa for all a G N n . And so, if L y is ^-continuous then sup{\y a \/w a : 
a G N n } < +00. ' □ 

Theorem 3.2. Let K C W 1 be as in IQj) and let y = a G N, 6e 

a given real sequence such that L y is £ w - continuous. Then y /ias a finite 
representing Borel measure fi on~K if and only if L y is nonnegative on Q(g). 
Equivalently if and only if: 
(3.1) 

L y (h 2 gj)>0 V/t € R[x], j = 0, . . . ,m; sup JM < M for some M. 

Proof. The necessity is clear. Indeed, if y has a representing measure fx on 
K then L y (h 2 gj) = J* K h 2 gjdfi > for all h G R[x] and all j = 0,1, ... ,m; 
and so L y (f) > for all / G Q(g). Moreover, ||L y ||^, < 00 because L y is 
^-continuous; hence (13. If) holds. 

Sufficiency. Suppose that L y is a non trivial £ w -continuous linear func- 
tional, nonnegative on Q(g), i.e., suppose that (]3.ip holds. In particular, 

L y {xf) < M(2k)\ for every k G N and every i = 1, ,n. Therefore, y 

satisfies Carleman's condition (12. 6p and since M&(y) ^ for all k G N, y 
has a representing finite Borel measure \x on M n , which in addition, is deter- 
minate. Next, using L y (h 2 gj) > for all h G M[x], and invoking Theorem 
12.21 one may conclude that gj > on supp fx, for every j = 1, . . . , m. Hence 
supp /x C K. □ 

Theorem 13.21 states that Q(g) solves the K- moment problem for £ w - 
continuous functionals. Of course, Theorem 13.21 is also true if one replaces 
the quadratic module Q(g) with the preordering P(g). 

The ^ w -closure of Q(g) and P(g). Observe that Psd(K) is £ w -closed. 
To see this, notice that with every every x G K is associated the Dirac 
measure <5 X , whose associated sequence y = (x a ), a G N™, is such that L y 
is ^-continuous. Indeed, let a := maxj|xj| so that |x a | < a'°', and let 
M := exp(a). 

M^ 1 \y a \ = exp(a)" 1 |x a | < exp(a)" 1 a |a| < a! < (2[|a|/2])! = w a , 

and so ||L y ||w < M, i..e., L y is £ w -continuous. Therefore, let (f n ) C Psd(K) 
be such that /|| w — ?■ as n — > 00. As L y is £ w -continuous one must have 
< lim^^oo Ly(f n ) = L y (f) = /(x). As x G K was arbitrary, / G Psd(K). 

Theorem 3.3. Let K C R n be as in / TO) and recall that Psd(K) := {/ G 
K[x] : / > on K}. Then cl w (P( 5 )) = cl w (Q( 5 )) = Psd(K). 

Proof. As Psd(K) is £ w -closed and Q(g) C Psd(K), cl w (Q( 5 )) C Psd(K), 
and so we only have to prove the reverse inclusion. Let / cl w (Q(<?)). 
Since Q(g) is a convex cone, by the Hahn-Banach separation theorem there 
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exists an £ w -continuous linear functional L-y that strictly separates f from 
cl w (Q(<?))- That is, there exists y £ N n such that L y is £ w -continuous, 
L y (f) < and L y (h) > for all h 6 cl w ((3(g)). By Theorem 13. 2\ such a y 
has a representing finite Borel measure fi on K, and so L y (f) = j K fd/j, < 
yields that necessarily /(xo) < for some xo G K. Hence Psd(K) C 
clw(Q(s , ))i which in turn yields the desired result. □ 

For instance, from Berg [3], the £i-closure of £[x] is Psd([— 1, l] n ). On the 
other hand, its ^ w -closure is now Psd(R n ), which is what we really want. 

It is worth mentioning that Theorem 13.21 and 13.31 can be extended to any 
preordering or quadratic module (finitely generated or not). The proof is 
the same. Any closed set K C R™ may be represented as the non-negativity 
set of such a quadratic module (taking generators of the form x i— > <?(x) : = 
Y17=i( x i ~ a i) ~ r2 > f° r suitable a £ R n and r > 0). However, for reasons 
that become obvious in the next section, the focus of the present paper is 
on the finitely generated case. 

4. Canonical £ w -projections onto Pa{g) and Qd(g) 

As we next see, the £ w - and £i-norm have the nice feature that one may 
find particular (canonical) projections onto various truncations of P(g) and 
Q(g) with a particularly simple expression. For £ w - and t\- projections to 
be well-defined we assume that K has a nonempty interior. 

We first provide an explicit form of canonical i\- and £ w -projection of 
any given polynomial / onto Pd{g) and Qd(g) respectively, and analyze their 
limit as d — > oo. Then we will consider the projections onto Pd(g) nl[x] s for 
fixed s, d G N, which, letting s —> oo, will permit to obtain the projection 
of / onto P(g) nR[x]i, and so to also characterize the sequential closure 

As in the previous section, for a polynomial in R[x]t we use indifferently 
the notation h for both the polynomial and its vector of coefficients h € R s ^. 
The context will make clear which one of the two is concerned. 

Let K C R n be as in (|2.2p and consider the following optimization prob- 
lem: 

(4.1) Pdw :=mf{||/-/ J || w : heP d (g)}. 

h 

That is, one searches for the best ^-approximation of / by an element h* 
of Pd{g), or equivalently, an ^-projection of / onto Pd(g)- In general, such 
a best approximation h* £ Pd{g) is not uniqu^3 but we provide a canonical 
solution with a very simple expression. Of course, and even though (|4.ip 
is well defined for an arbitrary / £ R[x], such a problem is of particular 
interest when / is nonnegative on K but not necessarily an element of P(g). 

lr The following example was kindly provided by an anonymous referee: Let n — 2, d = 1, 
x n> /(x) := — 2x\X2; and let C C R[x]2d be the cone of sums of squares of linear 
polynomials. Then any polynomial x t-¥ p\(x) := \(x\ — X2) 2 with A € [0,1], is an 
^w-projection (and an ^-projection) of / onto C and ||/ — Pa||w = 2. 
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Theorem 4.1. Let K C W 2 in 112. 2\) be with nonempty interior. Let / £ l[x] 
and let 2d > deg/. There is an £ w -projection of f onto Pd(g) which is a 

form: 



polynomial gj w G M[x]2d of the f 
(4.2) x ^ 5/ Pw (x) := /(x) +[X^ + Y / Y1 A ^ W 



,?;r 



2k 



(2*0! ; ' 



V i=l fc=i 

where the nonnegative vector A Pw € is an optimal solution of the 

semidefinite program: 

( n d n d 2k ^ 

^ik 

\ 'Mi t / / • / -r /mi -r / / 'V?. — 

A>0 



G ^G?) 



i=l k=l 



n d 



and Pdw = ||/ - ffflw = A Pw + ^ ^ 



i=l fc=l 



Proof. Consider / as an element of M[x]2<j by setting f a = whenever |a| > 
deg / (where |a| = X] a «)> an d rewrite (I4.ip as the semidefinite program: 



Pdv 



inf V] w a \ a 



(4.4) 



s.t. A a + /i Q >/ Q , VaGN^ 
K-K>-f a , VaeN^ 

m 



^ (Xj,B£> = 0, VqGN^ 

JC{l,...,m} 

A > 0; h € R[x] 2d ; Xj ^ 0, VJ C {1, . . . , m}. 
The dual semidefinite program of (|4.4p reads: 



sup ^2 fa( u 



(4.5) 



u,v>0,y 
S.t. 



aem 



u Q + v Q < w a Vq G 

Vq G N£ d) 



M rf (ffjy) h 0, VJC {!,..., m }, 



or, equivalently, 



sup -L y (f) 

s.t. M d (gjy) y 0, VJC{l,...,m} 
< ^a, Vq G N" d . 



(4.6) 



The semidefinite program (14.6P has an optimal solution y* because the fea- 
sible set is nonempty and compact. In addition, let y = (y a ) be the moment 

sequence of the finite Borel measure /J.(B) = J*Kn.B e ~" X " ^ or a ^ ^ G B, 
scaled so that \y a \ < w a for all a G N^. Then (y, u, v) with u = — min[y, 0] 
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and v = max[y, 0], is strictly feasible in (14. 5p . Indeed, as K has nonempty in- 
terior, then necessarily lVLj(gj y) >~ for all J C {1, ... , m}, and so Slater's 
condition^ holds for (|4.5p . Therefore, by a standard duality result in convex 
optimization, there is no duality gap between (|4.4p and (|4.5p (or (|4.6p ). and 
(j4.4H has an optimal solution (A*, (X*),^). Hence pd w = —L y *(f) for any 
optimal solution y* of (|4.6p . 

Next, recall that with J := 0, M d (# y) = M d (y). Moreover, M d (y) h 
implies M^(y) y for all k < d. By Lemma 1], Mjt(y) y implies 
that |y a | < max[L y (l), maxj L y (xf k )], for every a S N^,, and all k < d. 
Therefore, (|4.6p has the equivalent formulation 



(4.7) 



p d = -mi L y {f)) 
y 



s.t. m d {gjy) h 0, VJC {!,..., m} 



L y (l) < 1 



Indeed, any feasible solution of ()4.7p satisfies 

„2k 



L y {xf) < (2k)\, k = l,...,d; i = l,...,n. 



< max[L y (l),maxL y (xf )] < (2k)\ 



for every a with \a\ = 2k and 2k — 1, A: = 1, . . . , d. The dual of (|4.7|) is 
exactly the semidefinite program (|4.3|) . Again Slater's condition holds for 
(14. 7p and it has an optimal solution y*. Therefore (14. 3D also has an optimal 

solution A Pw G +1 with p dvi = A^ w + Yl=i Kk '» which is the 

desired result. □ 

The polynomial w € Pd{g) m (|4.2p is what we call the canonical -^w~ 
projection of / onto Pd(g)- 

Of course, all statements in Theorem 14.11 remain true if one replaces Pd(g) 
with Qd(g)- Moreover, if w a = 1 for all a (and so £ w is now the usual l\- 
norm) the polynomial w in (|4.2j) simplifies and is of the form: 

n 

(4.8) x ' y gf (x) := /(x) + (A p + £ Af x? d ), 

i=l 

for some nonnegative vector A p € R n+1 . If K = lR n then is the canonical 
•^-projection of / onto the cone of s.o.s. polynomials., as illustrated in the 
following example. 

Example 1. Let n = 2 and K = M?, in which case Pd(g) = Qd(g) = 
X[x]d for all d S N. Consider the Motzkin-like polynomial] x i— > /(x) = 



2 Slater's condition holds for the conic optimization problem P : inf x {c'x : Ax = 
b; x G K}, where K C R" is a convex cone and A £ R pxn , if there exists a feasible 
solution xo G int K. In this case, there is no duality gap between P and its dual P* : 
sup z {b'z : c — A'z G K*}. In addition, if the optimal value is bounded then P* has an 
optimal solution. 

Computation was made by running the GloptiPoly software [7] dedicated to solving 
the generalized problem of moments. 
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x\x\{x\ + x\ — 1) + 1/27 of degree 6, which is nonnegative but not a s.o.s., 
and with a global minimum /* = attained at 4 global minimizers x* = 
(±(l/3) 1 / 2 ,±(l/3) 1 / 2 ). The results are displayed in TableQ]for d = 3,4,5. 



d 


A* 


Pd 


3 


« 10- 


" 3 (5.445,5.367,5.367) 


« 1.610 -a 


4 


« 10" 


" 4 (2.4,9.36,9.36) 


« 2. 10~ 3 


5 


« 10" 


- 5 (0.04,4.34,4.34) 


« 8. 10~ 5 



Table 1. Best ^-approximation for the Motzkin po 



ynomial. 



4.1. Canonical £ w -projection onto fl l[x] 2 d and fl i[x] 2 d. 

We now consider £ w -projection of / onto Pt(g) H R[x]2<2 for given integers 
d,t G N, i.e., 

(4.9) pf w := inf { ||/ - 5 || w : a G P t ( 5 ) D M[x] 2d }. 

a 

In other words, we are interested in searching for a polynomial of degree 2d 
in Pt(g) which is the closest to / for the £ w -norm. For instance, if 2d = deg /, 
one wishes to find an ^-projection onto Pt(g) of same degree as /. One 
may also consider the analogue problem with the quadratic module, i.e., an 
^-projection onto Qt(g) H R[x] 2 rf. We also analyze the limit as t — > oo. 

Theorem 4.2. Let K C M n in \2. Sty be with nonempty interior and let 
d G N. Let f G R[x] and let 2t > max[2d, deg/]. There is an ^ w -projection 
of f onto Pt(g) nM[x]2d which is a polynomial g? w G M[x]2d of the fo 

/ n d od 

Pw x i 



form: 



(4.10) x ^ gf w (x) = /(x) + ^ + 5353 (,/,). y ' 

where the nonnegative vector A Pw G K nd+1 is an optimal solution of the 

semidefinite program: 

(4.11) 

( n d n d 2d 1 

- I i=l fc=l 8=1 fc=l V J 



and p dw = ||/ - a p Pw f || w = A Pw + £ £ A 



Pw 

i=l k=l 
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Proof. The proof is almost a verbatim copy of that of Theorem 14.11 with a 
slight difference. For instance, (14. 4p is now replaced with 



p?» := inf Y w a X c 
A,X,,/i ^ 



(4.13) 



s.t. X a + h a >f a , VaGN" d 
(4.12) K~ h a > -f a , Va G N% d 

h a - ( X J» B a>=0« Vtt€N^ 

JC{l,...,m} 

A > 0; /i G R[x] M ; Xj ^ 0, VJ C {1, . . . , m}. 
(where h a = whenever [a| > 2d) and the dual (|4.6p now reads 
sup -Ly(f) 

s.t. M ( ( ffJ y) >r 0, VJC{1 m} 

|Zfo| < Va G N" d . 

Again with exactly the same arguments, (|4.12j) has a feasible solution and 
(|4,13p has a strictly feasible solution y, and so Slater's condition holds for 
(I4.13p . which in turn implies that there is no duality gap between (|4.12|) and 
(|4.13p . and (I4.12j) has an optimal solution (A, (Xj), gf w )- However (and 
this is the only difference with the proof of Theorem 14. 1 p now one cannot 
guarantee any more that (|4.13p has an optimal solution because \y a \ < w a 
only for a G N 2(i and not for all a G N 2t . □ 

As before, we call gf w in (|4.10p the canonical ^-projection of / onto 
Pt{g) nM[x]2d. Of course, an analogue of Theorem 14.21 (with obvious ad hoc 
adjustments) holds for the canonical ^-projection <^? w onto Qt(g) nK[x] 2 d. 
Also and again, if w a = 1 for all a G N 2d , then the polynomial g^ w in (|4.10p 
simplifies to the form in (|4.8p . 

We next analyze the behavior of gj w as t — > oo to obtain the canonical 

£ w -projection of / onto P(g) nR[x] 2 rf. Recall that 



P(s)nR[x] 2d = \ (jPt(g) nM[x] 2d = \JP t d (g). 



t>0 I t>0 

Corollary 4.3. Let K C M n be as in \2. 2\i and with nonempty interior, 
f G R[x] 2d , and let g^{t) G P t (g) n M [x] 2( i be an optimal solution of U-9\) . 

Then there is an £ w -projection of f onto P(g) H M[x] 2d which is a polynomial 
gj G R[x] M of the form 

/ n d \ 

(4.14) x ^ gj (x) = /(x) + i\*o + E X *k *i d > 

V j=i fe=i / 

/or some nonnegative vector A* G M rfra+1 . In particular, z/K is compact and 
f > on K i/ien A* = and #J = /. 
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Proof. Let (pf*), t G N, be the sequence of optimal values of (14, 9p . which is 
nonnegative and monotone non increasing. Therefore it converges to some 
nonnegative value p dvf > 0. For every t G N, (14. 9p has an optimal solution 
9j- w (t) °f the form 

n d 2k 

(4.15) x i ^ 5/ Pw (t)(x) = /( x ) + A Pw (t) + ^^ Af fe w (t)|^, Vx € R», 

i=l fc=i 

with A Pw (t) > and Y,i,k X ik'( t ) = Pf ^ P^- Hence there is a sub " 
sequence (£;), j G N, and some nonnegative vector A* G R™ +1 such that 

3 j 



A Pw (tj) — > A* as j — > cxd. In addition, 



p dw = limpf = lim = A^^^Af*. 

1=1 k=l i=l k=l 

Hence gj w {tj) — > gj G K[x]2d as j — > oo where is as in (|4.14|) . And 
of course, as ||g pw (i) — gj\\ w —> 0, gj is in the closure P(g) nR[x]2d of 

f(j)ni[x] M . 

Next, suppose that there exists h G P(g) H R[x]2d with ||/ — h\\ w < \\f — 
gj\\ w . Then there exists a sequence (h t ) C R[x] 2( y;, t G N, with /i t G Pf{g) 
such that 1 1 /it — /i|| w — > as i — > 00. But then 

||/-/i|| w = \\f -ht + ht-h\\y, > \\f - ht\\w - \\ht - fejjw , Vt, 

>ll/-9f w Wllw ->° as t^oo 
and so taking limit as t — > 00 yields the contradiction 

||/-57!|w > [|/-/»||w > lim ||/-5f w Wllw = ||/-97llw 

J t— >oo J J 



Therefore ||/ - ^|| w = mm h {\\f - h\\ w : h G n R[x] 2d }. 

The last statement (when K is compact) follows from Schmudgen's Pos- 
itivstellensatz |15j which implies that if / is nonnegative on K then / + e G 
P{g) for every e > 0. □ 

Of course there is an analogue of Corollary 14.31 with Q(g) in lieu of P(g). 
The only change is concerned with the last statement where one needs that 
Q(g) is Archimedean. And also, if w a = 1 for every a G N n , then gj in 
(|4.14j) simplifies to the form (|4.8p . 

4.2. The sequential closures of P(g) and Q(g). Recall that for any 
convex cone A C R[x] 

(4.16) A x = {/ G R[x] : 3q G R[x] s.t. f + eq G A, Ve > }. 



We have seen that P(g) C P(g) x C P(<?), where ^4 denotes the closure of A 
for the finest locally convex topology. Therefore, it is of particular interest 
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to describe P(g)$, which the goal of this section. We know that 

(4.17) P{gf = |J P( 9 )nM[x] d , 

den 



and if for instance P(g) has an algebraic interior then P(g) = P(g)$. (See 
e.g. Kuhlmann and Marshall Prop. 1.4] and Cimpric et al. [5J Prop. 

1-3].) 

Notice that Q(g)* = \J deN Q(g) n R[x] d ([9]) and by [5, Prop. 1.3], we 
also have Q{g) = Q(g)^ if Q(g) is archimedean. 

Theorem 4.4. Assume that K C M. n in \2. 2\l has a nonempty interior and 
let f G R[x]. T/ien: 

(a) / G P(g)^ if and only if there is some d G N such that for every e > 0, 
the polynomial 

(4.18) x i-> f(x) + e(l + ^x 2 A ismP(g). 

(b) The same statement as (a) also holds with Q(g) instead of P(g). 

In other words, q G R[x] in J^.iffi ) for A = P(g) can be taken as x i— >■ 
{^ + YH=i x 'i d ) independently o/K. The dependence of q on f is through the 
power d only. 



Proof (a) From (j4.17j) / G P{g)$ if and only if / G P(g) n R[x] 2 d for some 
d G N. Next, by Corollary S3J / G P(p) nR[x] 2d if and only if the poly- 
nomial gj G R[x] 2c z defined in (|4.14p (and which simplifies to (|4.8p when 
ui Q = 1 for all a G N ra ) is identical to /. But then this implies that the 
polynomial gf w (t) G P t {g) n M[x] 2d in Kl5h is such that YJl=o A f W (*) ~> 
as i — >■ oo (recall that tt) Q = 1 for all a G N n ). Let \{t) := maxj[Af w (i)] so 
that X(t) — )■ as t — > oo, and the polynomial 



x ^ /(x) + A(t) ^ + E 
belongs to Pt(g) nR[x] 2 d because 



i=l 



f + X(t) 1 + *f = 5/ Pw (t) + A(i) - Aftt) + £ (A(t) - Af (t)) 
V / ' 25 ' ^ ' 



Therefore, for every e > 0, choosing t t such that A(i e ) < e ensures that the 
polynomial f+e(l+Y17=l x f d ) ls in -^< ! (fl , ) l ~ l ^[ x ]2d) which implies the desired 
result in (a). 

The proof of (b) is omitted as it follows similar arguments. Indeed, it 
was already mentioned that Theorem 14.21 and Corollary 14.31 have obvious 
analogues for the quadratic module Q(g). □ 
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Of course in f|4. 18j) one may replace x i w ^ n 1+X)£=i Y^k=i JWjl- 

In fact, and as pointed out by an anonymous referee, in (a) one may also 
replace 1 + Y17=i x i d with any point in the interior of Pd(g) in M[x] 2( i and 
in (b) with any point in the interior of Qd(g) in 1^ [ x ] 2d • The fact that 
1 + ^2i = ixf d is an interior point of Sfxj^ in M[x]2d (and hence also an 
interior point of Pd(g) and Qd(g) in M[x] 2 d) seems to be well-known. For 
instance, it can be deduced from [11] Theorem 3]. 

4.3. A Positivstellensatz for non compact K. As we know how to 
project with the £ w -norm, we are now able to obtain the following Posi- 
tivstellensatz on K. 

Corollary 4.5. Let K C W 1 in \2. 2\) be nonempty interior. Then f > on 
K if and only if for every e > there exists d G N such that 

(n d 2k \ 

Proof. The only if part: Recall that P(g) = Ud>o-^G?)' ano - from Theo- 
rem EES Psd(K) = cl w (P(#)). Let gf™(d) G R[x] 2d be the canonical £ w - 

projection of / onto P d (g) given in (g2J), where p dw = \o w +Ya=i Yfk=i ^ik* ■ 
As / G cl w (P(<?)), we necessarily have lim^oo p dvf = 0, because ||/ — 
gj w (d)\\ w — > 0. Hence given e > 0, let d be such that maxj ^ A^, w < e. 
Then 

/ n d ot, \ n d oju 

y. Il \ , „ z. n 



'^(2k)\ / ^ y ' v u ' ™ >(2k)\ 

1 fc=l K ' / i=l k=l V ; 



and so / + e (i + E « =i E J =i G P(5 ). 

The if part. Let </d G K[x] be the polynomial in (|4.19p . and let x € K be 
fixed, arbitrary. Then < Q'd(x) < /(x) + e ^T=i ex P Therefore, letting 
e -»■ yields /(x) > 0. □ 



5. Appendix 

Lemma 5.1. Let fx a finite Borel measure whose sequence of moments y = 
(y a ), a G N n , is such that for all i = 1, . . . ,n, and all k G N, L y (x? fc ) < 
(2fc!)M /or some M. Let f G R[x] 6e such that L y (xf /) > /or a// i = 
l,...,n, and aZZ i G N. Then the sequence = (<%), a G N n , where 
Za = L z f(x. a ) := L y (x a /) /or a// a G N™, satisfies Carleman's condition 

Proof. Let 1 < i < n be fixed arbitrary, and let 2s > deg/. Observe that 
whenever \a\ < k, |x| a < |xj| fc on the subset Wj := {x G W 1 \ [—1,1]" : 
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\xj\ = maxj \xi\}. And so, |/(x)| < ||/||i for all x € Wj. Hence, 

L zS (xf) = J /(x) xf d M (x) 

< / |/(x)|^ M (x) + *f + %(x) 
7[-l,l]« 

< M H/lli + MnH/ll! (2(fc + «))! < 2Mn||/||i (2(fe + a))!, 
and so we have 

L z/ (zfr 1/2 * > (2Mn||/|| 1 )" 1 / 2fc f((2(fe + a ))!)- 1 /2(*+«) N(A '"' ,) '" A ' 



> 



i(((2(A; + S ))!)- 1 /^) 



(k+s)/k 



1 / 1 
> - 



{k+s)/k 



2 V2(A; + s) 

where k > is sufficiently large so that (2Mn||/||i) _1 / 2fc > 1/2. Therefore, 

fc=i fc=fc v v ;/ 

where the last equality follows from (sp^T) ) (fc+s)/fe = ( 2(fc+s) )( 2(fc+ s ) ) s/fc and 
^ 2(k+s) — whenever fc is sufficiently large, say k > k\. Hence the 
sequence zf satisfies Carleman's condition f|2.6[) . □ 
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